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Some  Investigations  Relating  to 
the  Elastostatics  of  a Tapered  Tube 

by 

Barry  Bernstein 


1.  Introduction 

I 

The  problem  of  elastostatics  of  a tapered  tube  is  one  for 
which  one  searches  the  literature  in  vain.  The  problem  seems 
to  be  close  to  that  of  a cylindrical  tube,  but  this  appearance 
is  quite  deceptive.  Several  singularities  appear  in  the 
tapered  tube  which  do  not  in  the  cylindrical  tube.  Furthermore, 
in  a coordinate  system  appropriate  to  the  cylindrical  tube, 
separation  of  variables  is  possible.  Not  so  with  the  tapered 
tube.  In  this  report  we  shall  explore  some  approaches  to  the 
problem.  No  approach  that  we  shall  present  has  at  this  time 
shown  itself  to  be  the  clear  way  to  proceed.  However,  some 
of  them  may  have  some  promise. 

After  an  introduction  to  the  problem,  we  shall  discuss 
some  methods  found  in  the  literature.  Then  we  shall  present 
some  new  exploratory  results.  Although  we  cannot  be  conclusive 
at  this  time,  we  hope  that  we  have  opened  some  possibilities 
for  future  development. 

References  are  in  the  annotated  bibliography.  Section  8. 

Special  List  of  Symbols: 

Because  we  are  quoting  from  different  sources  which 
use  the  same  symbols  in  different  ways,  and  since  we  wish, 
with  only,  perhaps,  reasonable  modification  that  the  reader 
be  able  to  recognize  the  symbols  in  the  quoted  sources,  we 
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cannot  be  completely  consistent  in  using  a symbol  in  only  one 
way  in  this  report.  For  this  reason,  we  have  compiled  a list 
of  symbols  here  with  the  different  uses  of  the  same  symbol 
explained.  Symbols  are  listed  roughly  in  the  order  in  which 
they  appear  in  the  text,  except  that  all  listings  of  different 
uses  of  the  same  symbol  appear  together.  If  the  reader  will 
refer  to  this  list,  confusion  will  be  avoided. 


Symbol 

V 

V 

Cq,C2»C2,Cj 

V“ 

R,@  ,Z 
m 

ym 


Uses 

potential  function  (section  2) 
a solution  of  (28)  with  m=l  (section  7) 
constants  (section  2) 

Fourier  component  of  potential  (section  2) 
cylindrical  coordinates 
an  integer 

reduced  potential  (section  2) 


V2 

C =4  + i>^ 


c 


T,u,v 

u 


f 


defined  by  equation  (4) 

coordinates  for  the  tapered  tube  (section  2) 
(essentially  ^ and  of  section  2 are  the 
Niven  coordinates  j)  and  0 of  section  5) 

characteristic  coordinates  (sections  5,6) 

the  distance  of  the  intersection  of  the 
inner  and  outer  surface  of  the  tapered  tube 
from  the  axis 

separation  functions  (equation  6) 

some  arbitrary  function  (section  6) 

the  distance  from  the  singularity  in  Zak's 
coordinate  system  (section  3)  except  that 
Zak  takes  c = l,  this  is  the  same  as  the 
Niven  coordinate  r of  section  5 
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p 


n.r 

r 

r 

Fp(;^) 

X 

a 

f 

x,y,2 


r,e,^ 

e 

v,t 

¥ 

H 

fo*®o 

^0*^0 

F(a.P»'>f’,z) 


X,Y 


a Niven  coordinate  (sections  5,6,7)  in 
which  p = in  r 

Zak‘s  angular  coordinate  (section  3).  The  - 
same  as  the  Niven  coordinate  0 (sections  5,6,7 

Southwell  potentials  (section  3) 

gamma  function  (section  4) 

a curve  (section  6) 

a separation  function  (section  3) 

a constant  (section  3) 

the  taper  angle  (see  figure) 

a function  (equation  7) 

Cartesian  coordinates 

defined  by  equation  (23)  (section  6) 

same  as  V™,  but  in  Hein's  notation 

self  explanatory  - (equation  9) 

Niven  coordinates  (section  5) 
ti/2  - “S 

some  special  coordinates  (equation  15 
and  following  equation) 

some  constant  coefficients  (equation  17) 

defined  by  equation  (18) 

characteristic  function  (section  6) 

Niven  coordinates  of  a given  point  (section  6) 

characteristic  coordinates  of  a given  point 
(section  6) 

hypergeometric  function  (section  6) 

points  of  intersection  of  characteristics 
with  curve 

as  defined  in  the  equation  following  (25) 
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Q 


U 

fjjCz) 

P 

n 

h(n) 


as  defined  in  equation  (27) 

as  in  equation  (28)  - same  as 

functions  to  be  determined  (equation  29) 

an  index  (section  7) 

constants  to  be  determined  (section  7) 

an  integer  (equation  34) 

the  greatest  integer  in  n/2  (section  7) 

coefficients  to  be  determined  (equation  39) 


2.  The  Nature  of  the  Problem 

That  the  problem  of  elasto statics  hangs  on  the  study  of 
Laplace's  equation  is  well  known.  A review  of  solutions  of 
such  problems  in  terms  of  potential  functions  is  given  by 
Green  and  Zema,  section  (5,6)  [13].  If  one  could  handle 
Laplace's  equation  for  the  tapered  tube,  then  elastostatic 
problems  would  be  accessible. 

The  first  effort,'  then,  that  seems  reasonable  is  to  see 
if  separation  of  variables  is  possible.  We  turn,  then,  to  the 
work  by  Snow  [l].  We  consider  here  chapter  IX,  p.  228  of  this 
work. 


In  cylindrical  coordinates  R,  Z,  ® , we  have  for  Laplace's 
equation  for  a potential  V 


^ 1 av  ^ a£v_  _ a 

"T  - + A Q — VJ  • 

R aR  R'^  3®'=^ 


(1) 


From  equation  (1),  we  may  immediately  separate  out  the 
angular  coordinate  ® by  writing  V in  a Fourier  series  in 
Indeed,  Slow  writes 
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V = Cq  + C^z  + (C2  + C^z)  inR 

+ T v“(R,Z)co5m  (®  - 
m=0  “ 


where  Cq,  C^,  C2»  and  are  constants.  The  coefficient 

satisfies 


dV  1 ^ _ m2v®  ^ 

^ ^ ^ 3r  ”r^ 


(2) 


Or,  putting 

1 

V“  = R ^ U“ 


in  equation  (2)  one  obtains  an  equation  for  the  reduced 
potential  u“,  namely 


(3) 


where 

W 

aR^  az^ 

Now  if  one  looks  at  the  R-Z  plane  one  sees  that  the  trace 
on  this  plane  of  a tapered  region  can  be  represented  as  a wedge, 
one  side  of- which  is  parallel  to  the  Z-axis  at  some  distance, 
say,  c,  from  the  axis,  and  crossing  the  other  side  at,  say, 

the  Z-axis  at  an  angle  a (see : figure) , which  we  call  the  , 

\ 

taper  angle.  A conformal  mapping,  then,  from  the  Z + iR  plane 
into  the  plane  of  T = ^ + i7^  given  by 

C=  log(Z  + R - iC) 
gives  for  (3) 
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(5) 


+ 


0. 


Now  equation  (5)  is  an  equation  in  4 . , which  are 

natural  coordinates  for  the  tapered  tube.  Indeed,  in  this 
coordinate  system,  the  surfaces  of  the  tube  become  C = 0 and 
^ = a.  One  may  say  more;  This  is  essentially  the  only  ortho 
gonal  coordinate  system  in  which  the  surfaces  of  the  tube 
become  coordinate  surfaces  for  any  value  of  a. 

We  now  ask  the  question  whether  or  not  one  may  find 
solutions  of  the  form 


v“  = T(4  ,)^)u{^)v(rj)  (6) 

where  T-  is  to  be  found.  The  answer  is  given  by  Snow  (pp.  252- 
253).  It  appears  that  the  answer  is  no,  since 


e^^ 

' (a  + e^  sini^ 


6e^^  cos^  (l  + e^  cosr^ ) 
(a  + e^  sinTj 


^ 0 


which,  by  application  of  Snow’s  result  to  our’  equation  implies 
that  separation  of  variables,  even  to  within  a known  factor  T, 
is  not  possible  for  the  coordinate  system  ( 4 » >^  )• 

Techniques  of  separation’  of  variables,  with  all  their 
ramifications,  then  fail.  Other  techniques  must  then  be  sought. 
And  a look  at  some  of  these  is  then  our  task. 


3.  The  Work  of  Zak 

Here  we  shall  discuss  a technique  used  by  Zak  for  solving 
a problem  of  a cylinder  with  stress  singularities.  The  method 
happens  to  involve  the  Southwell  potentials  [l4],  but  ’the 
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essential  feature  of  it  is  the  method  of  obtaining  a sequence 
of  functions  which  approach  a solution. 

Zak  used  a coordinate  system  which  is  essentially  that 
developed  in  the  previous  section  and,  indeed,  is  equivalent 
to  Niven’s  coordinates  (section  5)*  If  we  replace  ^ by  the 
letter  and  ^ by  the  letter  we  shall  have  the  coordinate 
system  which  he  uses.  In  this  section  we  shall  adhere  to  Zak’s 
notation.  Hov/ever,  the  same  letter  ^ will  be  used  differently 
elsewhere  in  this  report,  so  caution  on  the  part  of  the  reader 
is  urged.  Please  refer  to  the  list  of  sjmbols.  If  referred 
to  the  tapered  tube,  Zak's  coordinates  are  the  Niven  coordinates 
normalized  so  that  c = 1. 

The  Southwell  potentials  as  modified  by  Zak  satisfy  the 
equations 


^ 1 dn  ^ 

^ dR  a?” 

= 0 

1 ar  a^r 

^ a^n 

aR^  ' ^ aR  ' az^ 

az^ 

in  cylindrical  coordinates.  After  expressing  these  equations 
in  terms  of  his  and  Zak  seeks  a solution  for,  say,  /I  in 
the  form 

p=0  I P 

and  obtains  a sequence  of  equations 

Fo"  + = 0 


n 


-cos/zf  j 

so  that  each  f;mction  F depends  on  •ttie  previous  ones.  A 

It 

similar  technique  is  applied  to  P. 

It  is  not  difficult  to  see  that  Zak's  technique  coiild 
readily  be  applied  to  the  tapered  tube  problem:  The  proper 
coordinate  form  and  the  technique  are  already  developed. 

At  the  time  of  writing  of  this  report,  we  feel  that  Zak's 
method  may  be  the  most  promising  where  it  can  be  applied.  It 
appears  to  have  two  disadvantages.  Zak  expands  a term  as 

■ mfo'  f 

which  has  as  its  domain  of  convergence  a region  near  p = 0. 
(This  region  was  of  interest  for  the  study  of  a singularity 
at  j)  = 0. ) 

For  the  tapered  tube  it  may  be  of  Interest  if  one  limits 
oneself  to  regions  where  1^  sin;^|  < 1,  but  this  means  that 
the  radial  length  allowed  is  limited  by  the  angle  of  taper. 

For  example,  for  a taper  angle  of  2®,  the  expansion  is 
valid  for  ^ up  to  about  28  (distance  from  singularity  about 
28  times  the  quantity  C in  Niven  coordinates)  and  convergence 
would  probably  be  slow  if  j)  were  near  28. 

Although  we  do  not  see  how  to  do  it  at  present,  it  may 
be’  possible  to  apply  Zak's  technique  to  a far  away  region. 
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But  the  trouble  at  the  moment  is  that  as  one  goes  toward 
larger  ^ there  are  points  closer  and  closer  to  the  surface 
^ = 0 at  which  (1  - p sin^z^)”^  becomes  infinite. 

The  second  disadvantage  which  may  be  minor  is  that  one 
does  not  deal  with  a sequence  of  exact  solutions.  However, 
this  would  not  necessarily  impair  its  usefulness  where  con- 
vergence is  rapid  enough.  Nevertheless,  in  the  broad  study  of 
the  question,  ‘a  search  for  exact  solutions  should  be  made.  If 
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problem,  they  might  or  might  not  provide  a better  method  than 
that  of  Zak  in  some  given  situation.  In  sections  6 and  7 we 
report  on  a search  for  such  solutions. 


4.  Integral  Equations 

The  method  of  integral  equations  rests  on  the  representa- 
tion of  the  solution  of  Laplace's  equation  as  an  integral.  A 
review  of  such  integral  representations  is  given  by  Temple  [7]> 
who  contends  that  the  culmination  of  this  work  is  in  Whittaker's 
result  that  potential  functions  which  are  regular  near  the 


^Jx  iiovc  j.  ua 


!2ti 

f (z  + i X cos  p + i y sin  p,p)dp. 

0 


(7) 


Basically  the  method  of  integral  equations  consists  of  setting 
up  equations  for  the  unknown  function  in  an  integral  expression 
such  as  (7).  ‘.These  equations  are  based  on  the  bomdary 


conditions. 


A review  of  the  use  of  the  method  of  integral  equations 
is  given  by  Heins  [2],  who  makes  use  of  the  Poisson  integral 
representation:  For  a function  Sj^(R,Z)  satisfying 


d^i. 


n 


dZ" 


(8) 


which  is  the  equation  ' satisfied  by  a Fourier  component 

or  of 
n n 

*. 

$(R,/zf,Z)  = iQ(R,Z)  + f S^(R;Z)  cosn)!^ 

n=l 


00  _ 

+ Z sinn^ 

n=l 


(9) 


of  an  harmonic  function  i,  one  obtains 

.Tt 


U„(R,Z)  = — f U (0,z+ ir  cosY  )sin^nYd"t 

r(4)  r(n+  f)  Jq  ^ 


where 


= i. 

n n 


(10) 

(11) 


Now  the  validity  of  the  Poisson  Integral  Representation 
(10)  hangs  upon  the  regularity  of  the  solution  on  the  Z axis. 
Indeed  the  assumption  of  such  regularity  is  stated  explicitly 
by  Heins  (p.  789)  and  the  problems  solved  (e.g.  a charged 
disc,  or  a lens,  with  axis  along  the  Z-axis)  do  not  violate 

Note  Equation  (8)  is  the  same  as  equation  (2)  using  Heins' 
notation  instead  of  Snow's. 
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this  condition.  Other  work  which  we  have  found  so  far  [9,  10, 
11,  12]  does  not  seem  to  violate  this  condition. 

¥e  must  caution  that  we  have  not  at  the  time  of  writing 
of  this  report  fully  digested  the  question  of  whether  regularity 
on  the  Z-axis  is  absolutely  crucial  to  whether  or  not  the 
problem  of  the  tapered  tube  is  amenable  to  a Poisson  Integral 
type  analysis.  However,  the  Z-axis  is  outside  the  domain  of 
required  validity  of  solutions  to  such  problems.  So  there  is 
no  reason  to  expect  iiiat  the  Poisson  Integral  will  give  the 

I 

answer.  On  the  other  hand,  neither  can  one  assert  at  this 
point  that  it  will  not  figure  in  a method  of  solving  the  tapered 
tube  problem.  Indeed,  perhaps  we  need  a solution  valid  outside 
the  inner  surface  as  well  as  a solution  valid  inside  the  outer 
surface  of  the  tapered  tube,  so  that  their  region  of  common 
validity  will  be  as  desired. 

Another  method  which  we  feel  needs  further  exploration 
is  that  of  Snow  [l],  Chapter  IX.  Again,  we  feel  at  the  time 
of  writing  of  this  report  that  we  have  not  yet  seen  through 
the  method  well  enough  to  be  certain  that  it  will  apply  in 
whole  or  in  part  to  the  tapered  tube.  The  difficulty  at  the 
axis  arises  in  trying  to  map  the  R-S  plane  into  the  wedge- 
region  which  is  that  of  the  tapered  tube  on  the  R-Z  plane 
without  getting  into  the  same  type  of  difficulties  with  the 
axis.  However,  for  reasons  similar  to  those  stated  in  con- 
nection with  the  Poisson  Integral,  we  feel  that  the  matter  is 
not  at  all  settled  at  this  time  and  that  we  should,  indeed, 
like  to  consider  it  further. 
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Nevertheless,  In  order  to  seek  integral  equation  solutions 
appropriate  to  the  tapered  tube,  it  would  be  nice  to  have  an 
integral  representation  which  is  tailored  to  hold  in  the  proper 
region.  To  this  end,  we  have  carried  out  an  investigation 
based  on  the  theory  of  characteristics.  It  may  seem  odd  to 
do  this  today,  but  in  nineteenth  century  work,  the  relation  of 
the  wave  equation  to  Laplace's  equation  through  the  use  of 
complex  characteristics  was  well  accepted.  We  shall  present 
these  results  as  soon  as  we  have  discussed  the  work  of  Niven. 

5.  The  Work  of  Niven 

A coordinate  system  appropriate  to  the  tapered  tube  was 
treated  by  Niven  [4].  Indeed,  he  defines  a coordinate  system 
r,  ©,  ^.by 

X = (C  +rcose)cos/^ 
y = (C  + r cos9)sin^ 
z = r sin  0 

(where  the  tilda  is  our  notation). 

We  find  it  more  convenient  to  deal  with  the  complement 
of  0.  Thus,  we  shall  interchange  sin0  and  cos0  to  write 
X = (C  + r sin0)cos^ 
y = (C  + r sin0)sin)^  . 

2 = r COS0 

and,  since  these  differ  so  trivially  from  Niven's  coordinates, 
we  shall  call  these  Niven  coordinates  also.  It  is  clear  then 
that  in  the  Niven  coordinates  the  surfaces  of  the  tapered  tube 
with  taper  angle  a are  simply  0 = 0 and  0 = a. 
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Niven  then  presented  a form  of  Laplace/ s equation  in 
these  coordinates,  namely 

— r (c  + r cos9)  ^ - -4;  (c  + r cosG)  ^ 

8r  Sr  89  39 


3^V 


= 0 . 


c + r cos9 

It  is  clear  that  by  making  the  substitution  9 = n/2  - 9,  one 
gets  instead 

— r (c  + r sin©)  — + — — (c  + r sin©)  — 

3r  dr  3©  9© 


B^V 


= 0. 


c + r sin©  bf^ 

Niven  then  writes  this  equation  in  terms  of  p = ^n  r for 
the  case  where  V is  independent  of  We  shall  write  instead 
the  equation  in  the  general  case,  namely 

— (c  + ePsin©)  ~ — (c  + e^^sin©)  — 

Bp  Bp  B9  39 


.2p 


B^V 


= 0. 


c + ePsin©  3^^ 

These  Niven  coordinates  have,  in  fact,  been  met  before 
in  this  report  and  indeed  p and  0 are  respectively  the  ^ and 
of  section  2.  In  this  section,  however,  we  shall  use  ^ 
and  for  other  coordinates. 

Niven  writes  Laplace's  equation  for  the  case  where  V is 
independent  of  ^ in  terms  of  (characteristic)  coordinates 


. p+i© 

s = c + eJ  , 


. ef 


-i© 


(13) 
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as 


2(<*  + T^) 


0. 


(14) 


To  him  this  is  Just  a step  in  a process  which  in  fact 
does  not  seem  to  lead  us  very  far.  Indeed  the  next  step  is 
to  put 


4 + 

= t 


(15a) 

(15b) 


and  arrive  at 


e2v  £v  _ 0 
dt^ 


(15) 


which  does  indeed  admit  separation  of  variables  solutions. 
However,  the  simple  observation  that 


V = log2R,  t = 2i2 

shows  that  indeed  we  are  merely  back  in  a slightly  unfamiliar 
form  to  cylindrical  coordinates  R and  Z and  that  the  separation 
of  variables  will  yield  the  familiar  Bessel  functions. 

The  rest  of  Niven's  results  are  of  mathematical  interest, 
but  do  not  help  us  with  our  tapered  tube  problem.  He  does 
indeed  obtain  some  closed  form  solutions,  but  they  do  not  have 
the  proper  equipotential  surfaces  appropriate  to  the  tapered 
tube  problem.  Indeed,  his  method  is  to  seek  solutions  of  the 
form 

2 A e^'^^cosn©  (16) 

n=0 

where  m depends  on  n in  some  way,  and  then  cleverly  to  choose 
some  such  dependence  which  allows  the  series  to  yield  a solution. 
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The  method  is  interesting,  but  it  probably  lacks  the  generality 
which  we  need  in  order  to  obtain  enough  solutions  that  we  can 
handle  the  tapered  tube  problem.  Indeed,  in  section  7 below 
we  shall  present  solutions  containing  types  of  terms  not  fovind 
in  (16). 

In  closing,  we  point  out  that  these  type  of  coordinates 
were  used  by  Riemann  [7]  to  solve  a problem  of  an  anchor  ring. 
Although  Riemann *s  paper  is  very  pretty,  again  it  contains  no 
hint  how  to  seek  solutions  with  equipotential  sxirfaces  appro- 
priate to  the  tapered  tube. 

6.  An  Integral  Representation  of  the  Solution  to  Laplace's 
Equation  in  Niven’s  Coordinates 

Consider,  now,  the  form  of  the  equation  for  the  potential 

V^^^(R,Z)  in  polar  coordinates  (R,  <9)  as  given  by  equation  (2). 

We  shall  obtain  formally  an  integral  representation  of  the 

solution  to  this  equation  in  the  Niven  coordinates  (j),0),  where 

R = c + e^sin©  (17a) 

Z = ef^cos©.  (I7b) 

To  begin  with,  let  us  put 

= r“w,  (18) 

t 

(where,  of  course,  W will  also  depend  on  m,  but  we  find  it 
notationally  simpler  not  to  write  it  explicitly).  We  obtain, 
then  for  (2) 

v2w  . ^ . 0.  (19) 

aR 

Now  let  us  change  to  coordinates  4 , ^ » where 
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4 = R + iZ  = C + ieP'^® 

(20a) 

T|=R-iZ  = C-  ieP"^^®. 

(20b) 

(Note  that  (20a)  and  (20b)  follow  from  (17)  and  (13). ) 
We  then  obtain  for  (19) 

. 2in+l  1 . dW 

+ — i5 — T I TT  + 


)■ 


0. 


(20) 


Now  (20)  'is  in  characteristic  form,  albeit  the  char- 
acteristics are  complex,  either  in  the  R,  Z coordinates  or  the 
p,  6 coordinates.  It  is  clear,  then,  that  the  characteristics 
are  given  by 

p + i©  = constant  (21a) 

p - i©  = constant.  (21b) 


and 


Now  equation  (20)  is  of  the  form  given  in  Sommerfeld's 
book  [15],  section  11,  in  relation  to  a hydrodynamic  example 

+ VvT^  D -1  Vi  t Q ti+’t  /Mn  +-/'\ 

UX  l^y  i UlLXV/VI  U V.*.  I^.J.Xh'XX  - 


adjoint  equation 


2nH-i  a r 1 pH  ^ aH\i 

bl  L4+71 


= 0 


for  each  point  ^©(po’^O^  such  that  H = 1 at  P and 
SH  2m+l 


2(^  + 1^) 


H = 0 


and 


b>l  2(4  + )^) 
on  the  characteristics 
which  pass  through  pQ,  ©q. 


H = 0 


const  and  4 = const  respectively 
Riemann  solved  this  problem.  The 
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solution  for  our  case  is 


where  . 

(4  - - ’Iq) 

P(a,p,y,2)  is  the  hypergeometric  function 


(22) 


(23) 


F(a,p,T,z) 


1 + 


^ 2 

T T7 


a(a-Hl)p(B+l) 
YCr  +1) 


and 


4o  - c 
^0  = ^ 


. >-i®0 

+ 1 e 

- 1 e 


(24a) 

(2Ab) 


To  obtain  an  integral  form  of  the  solution  of  (20),  then 
we  apply  the  relation  (7)  p. (54)  of  Sommerfeld  [15],  which  we 
write  -in  the  form 

= j (X<i4  - Yd7|)  + ^(WH)p^  + ^(WH)p^  (25) 


where  T is  a portion  of  a curve  on  which  data  are  given, 

and  P2  are  the  intersections  of  the  curve  with  the  respective 
characteristics  through  (Pq>6q)»  "the  integration  is  taken  to  be 
along  r from  P^^  to  P2*  (Here  we  must  be  careful  that  P^^  be 
the  point  where  TJ  = intersects  P and  P2  is  the  point 
where  4 = 4q  intersects  P. ) Also 
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2m+l 


2 \ di|  / 

Y WM) 

2 \ / 


2(4  +i|) 

2m+l 
2(4 + TJ) 


HW 


HW. 


We  wish  now  to  evaluate  the  right  hand  side  of  (25), 
where  we  shall  take  for  P the  curve  6 = 0,  which  corresponds 
to  that  surface  of  the  tapered  tube  which  is  parallel  to  the 
axis.  We  caldulate  first  for  any  function  u 


Bx  2 L 

du  _ e^P"^^  r i ^ + Bu~| 

9y  2 L 

d4  = i e?”^®(dp  - i d0) 

dT^  = -i  eP''’^®(dp  +id0). 

Then  we  obtain  from  (25)  and  the  relations  following  it 

rPo^i  ®0 

V(pO,©o)  = i 1 Q(y,0)do  + 

Po"^  ®0 

|(HW)(^Q-ieQ,0)  + |(HW)(pQ  + i 0Q,O)  (26) 

where 

Q(p,e)  = ? [^  H]  ^27) 


Now  (26)  and  (27)  give  an  integral  representation  of 
the  solution.  The  limits  pQ  - i 6q  and  pQ  + i ©q  are  of  course 
the  (complex)  values  of  p at  which  the  characteristics  (21a) 
and  (21b)  through  (po»®o^  intersect  6 = 0. 
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We  have  been  proceeding  formally,  but  it  is  clear  that 
we  must  assume  analytic  boundary  data,  for  otherwise  the 
integrand  will  be  path  dependent. 

We  shall  argue  that  the  Integrand  Q(j),0)  is  real  for  real 
values  of  Assume  this  for  the  moment. 

If  we  make  the  change  of  variable 

p » po  + 

in  the  integrand,  we  get  for  the  integral  term 


By  putting  - ^ for  ^ in  the  first  integral  on  the  right 
hand  side,  we  get  for  the  integral  term 


- j [Q(p  + * Q(f  - 1/,0)W, 


which,  under  our  assumptions,  becomes 
f®0 

-2  ReQ(pQ  + l^,0)d;2^. 

We  obtain, then 

W(p,e)  = 

.0 


Re 


+ H(p  + ie,o)w(p  + ie,o) 
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where  the  Integrand  is  evaluated  at  (|)  + i)z^,0).  This  is  a form 
of  the  integral  representation  of  the  solution. 

Now  let  us , look  at  the  assumptions  on  Q.  First  of  all , 

t 

in  order  to  be  sure  that  our  operations  are  legitimate,  we 
must  assign  analytic  bomdary  data  on  W(^,0)  and  dW/d0  at 
(y»,0).  These  must  be  assigned  so  that  they  are  real  when  p 
is  real.  Next  look  at  H.  There  is  no  problem  of  analyticity 
of  H.  However > we  shall  check  that  H is  real  when  6 and  ^ are 
real,  so  that  it  will  also  follow  that  dH/d0  will  be  real,  at 
real  (j),0). 

¥e  have  from  (20)  and  (24)  that 
4 + I'l  = 2(c  + efsin  0) 

^0  ■*’lo  “ 

Now  the  last  equation  is  the  product  of  an  expression  with 
that  of  its  complex  conjugate  and  is  thus  real.  From  (23),  it 
then  follows  that  z is  real  for  real  (p,0),  and  hence,  from 
(22)  and  (23)  we  see  that  H is  real  for  real  (p,0). 

This  establishes,  at  least  formally,  the  integral  repre- 
sentation. 

7.  Some  Polynomial  and  Series  Solutions 

Although  it  is  not  possible  to  separate  variables  in  the 
Niven  coordinates,  or  in  an  equivalent  coordinate  system,  it 
is  possible  to  obtain  some  exact  solutions  of  more  generality 
than  those  given  by  Niven.  To  this  end  let  us  start  with 
equation  (3)  for  the  reduced  potential  in  cylindrical 
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coordinates.  We  drop  the  superscript  m for  notatlonal  con- 
venience. We  have,  then 


R^V^U+  (|  - in^)U  = 0.  (28) 

We  substitute  for  U the  series 

U = I fAZ)  (29) 

k=0  ^ 

where  the  index  p and  the  functions  fjj(Z)  are  to  be  determined. 
Substitution  of  (28)  into  (29)  yields 

[p(p-l) +^-m2]fQ(Z)  + [(p+l)p+|-m2]f^(Z)R 

+ ^^{[(p+k)(p+k-l) +J-m^]fj^(Z) +fj^'_2(Z)]  R^  = 0.  (30) 

Setting  the  coefficient  of  fQ(Z)  to  zero,  we  get  the 
indicial  equation 

n 2 O 

(p  - ^)  - = 0 

or 

p = ^ +in. 

We  choose  the  larger  root  p = ^ + m of  the  indicial  equation. 
(As  with  Bessel  functions,  the  smaller  root  will  give  nothing 
new. ) We  obtain  by  substitution  into  (30) 

(2m+2)fj^(2)  = 0 (31) 

k(2m+k)fi^(Z)  + f»  (Z)  = 0.  (32) 

^ k-2 

Equations  (31)  and  (32)  taken  together  tell  us  that 
fj^(Z)  = 0 for  odd  k 
and  that  for  even  k,  k = 2j,  we  have 
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(33) 


2"^0(0+m) 


This  gives  us  a method  of  obtaining  solutions,  since  we 
can  now  pick  fQ{<S>)  as  we  wish  and  make  use  of  the  recursion 
relation  (32).  We  shall  explore  some  of  the  consequences  of 
this  procedure. 

We  see,  in  particular,  from  (33)  that  if  we  choose  fQ  to 
be  a pol3niomial  in  Z,  then  the  number  of  terms  in  the  expansion 
(29)  will  be  finite.  The  solution  will  be  a square  root  of  R 
times  a polynomial  in  R and  Z.  We  shall  Investigate  the  case 

,fQ(Z)  = Z^,  n = 0,1,2,....  (34) 

From  these  one  can,  in  fact,  by  addition  and  multiplication 
by  constants,  arrive  at  the  results  for  any  choice  of  a poly- 
nomial for  f^. 

Now  (34)  and  (33)  tell  us  that 

^23  * . <55) 

for  some  Cj  to  be  determined.  If  we  put  (35)  into  (35)  we 
obtain 


C,  = . ^ 


The  solution  to  the  recursion  relation  (36)  is 


(36) 


C,  . (-1)3  -JaiJE O- 

^ 2^3  O- 

where  j = 0,1,...  ,n/2  if  n is  even  and  j = 0,1,. 
n is  odd. 


(37) 

. ,(n-l)/2  if 
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If  we  take  c^nlm!  = 1,  we  obtain  then  a solution 

„ R2d2n-2d  (Jg) 

where  h(n)  = n/2  if  n is  even  and  h(n)  = (n-l)/2  if  n is  odd. 

Equation  (38)  gives,  then,  a square  root  of  R times  a 
polynomial  in  R and  Z.  Now  returning  to  Niven  coordinates, 
we  obtain 

U=  (c+r  sine)”^^^(r  cose)“^^  (^sece+tanG)^^. 

j=0  2 ^(n-2o).'o'.  (n+j)! 

Since  there  are  finitely  many  terms,  there  is  no  problem 
of  convergence. 

Another  set  of  exact  solutions  can  be  generated  by  writing 
f^CZ)  = z-^ 

where  n is  an  integer.'  We  obtain,  then,  by  the  same  process 

0=0  2^^  d!  ( j+m)i 

According  to  -the  ratio  test,  this  series  will  converge 
for  1^1  < 1. 

In  the  Niven  coordinates,  then  we  obtain 

U = r""(o+rsine)^  J (£secettane)^J 

J=0  2^0 

which  converges  when 

c+r  sec©  < rcos©  (1st  quadrant) 
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cos©  - sin©  ’ 


Now  let  us  consider  a taper  angle  a < ^,  so  that  we  are 
concerned  with,  say, 

0 $ © ^ a. 

Then  cos©-  sin©  is  decreasing  with  ©,  since 


;( cos© -sin©)  = -sin© -cos©. 


Thus 


cos©  - sin©  cosa  - sina 


, 0 ^ © ^ a. 


Therefore  convergence  in  the  wedge  region  will  occur  for 


cosa  - sina 

which  means  that  the  solution  will  be  valid  away  from  the 
singularity  where  the  outer  and  inner  surface  of  the  tapered 
tube  meet.  This  is  the  sort  of  region  that  figures  in  a gun 


tube. 


For  the  case  of  m = 0,  we  have  also  succeeded  in  finding 


an  additional  solution  involving  a logarithm  and  a polynomial. 
We  seek  solutions  to  (28)  in  the  form 

00  i+2j  „ 

U = V in  R + J b.  (39) 

0=0  ^ 

where  V satisfies  (28)  with  m = 0.  We  obtain 


0 = RV%  + g = [R^72u,+  g]in  R 


+ 2r|^-V+  I (4a^-i) 

Sr  j=o  ^ ^ 


29 


We  take  for  V the  solution  (38)  which  we  write  as 


V = 


where  c.  is  given  by  (37)*  After  some  straightforward  manipu- 
lation  we  obtain  the  recursion  relation 


^[(n-2j+2)(n-2a+l)b^^^  + 4jCj],  j = 1,2,... 


(40) 


We  note  that  since  there  are  only  a finite  number  of  non-zero 
Cj,  and  since  the  coefficient  of  (40)  becomes  zero  for 

j = (n+2)/2  or  j = (n+l)/2  (depending  on  whether  n is  even  or 
odd)  there  are  only  a finite  number  of  non- zero  b.  and  so  both 
V and  the  summation  term  in  (39)  become  polynomials. 

Hiese  solutions,  then,  in  Niven  coordinates  become 


(n+i)/.^ ^ — 

' K ^ + Sin© 


h(n+2) 

L ^0 


[^sec0+tan0]^^[c  .in(c+r  sin©)  + b.] 


Some  thought  must  be  given  to  the  best  way  to  assemble 
these  solutions  in  order  to  satisfy  given  boundary  conditions. 
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This  work  seems  very  rich  and  may  well  hold  the  clue  to 
a procedure  for  attacking  the  tapered  tube  problem.  It  does 
establish  criteria  for  whether  one  may  separate  variables  for 
Laplace's  equation  after  transforming  the  polar  coordinates 
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This  is  a review  of  the  use  of  the  Poisson  Integral 
formulation  to  solve  problems  of  potential  theory  by  the  means 
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solution  is  represented  by  an  integral  whose  integrand  involves 
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such  as  the  Helmholtz  representation,  for  example,  is  then 
used  to  set  up  a relation,  actually  an  integral  equation,  with 
which  to  determine  the  unknown  integrand  in  the  Poisson  integral. 
When  this  relation  can  be  solved,  then  one  has  the  solution 
to  the  problem  represented  as  an  integral.  Problems  given  as 
examples  include  the  electrostatic  fields  about  discs  and 
circular  lenses  with  axes  on  the  Z axis. 

It  is  not  clear  at  the  time  of  writing  of  this  report 
whether  or  not  the  method  can  be  effectively  used  for  the 
tapered  tube  since  1)  the  method  is  basically  designed  for 
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functions  which  are  regular  on  the  Z axis  (or  at  least  a 
portion  of  it)  and  2)  its  application  seems  to  depend  on  its 
proper  use  in  an  appropriate  coordinate  system.  Whether  and 
how  to  use  it  in  whole  or  in  part  for  Idle  tapered  tube  problem 
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with  Stress  Singularities",  J.  Applied  Mech.  39  ( E ) 
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the  Niven  coordinates.  The  method  does  not  use  a sequence  of 
exact  solutions,  but  it  does  obtain  a series  which  converges  to 
an  exact  solution  where  it  converges.  As  the  method  is  applied, 
it  is  limited  by  a singular  curve  (psin/=  1),  a function  which 
occurs  in  the  equations  (1  -p  sin;!()"^.  As  the  distance  p from 
p = 0 (which  would  be  the  intersection  of  the  inner  and  outer 
surfaces  in  the  tapered  tube  problem)  gets  larger,  the  location 
of  the  singular  curve  gets  closer  and  closer  to  where  sin^  = 0, 
which  limits  the  region  in  which  Zak's  method  could  be  used,  at 
least  without  some  modification.  We  calculated  that  for  a taper 
angle  of  2*  one  is  limited  to  a distance  less  than  about  28  times 
the  distance  from  the  axis  to  the  intersection  of  the  inner  and 
outer  surfaces  of  the  tapered  tube. 
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9.  Conclusion 

The  problem  of  the  best  way  to  proceed  to  obtain  analytical 
solutions  of  the  tapered  tube  problem  remains  unresolved  at  the 
moment.  We  feel  that  some  progress'  has  been  made  in  developing 
exact  solutions  in  the  appropriate  coordinate  system.  Perhaps 
a method  of  assembling  these  solutions  to  fit  boundary  data 
would  be  fruitful.  Perhaps  a formulation  in  integral  equation 
form,  would  be  the  way  to  proceed.  Again  we  feel  that  some 
progress  has  been  made  by  our  integral  representation  of 
section  6.  And  of  course  Zak‘s  basic  idea  could  be  developed 
further.  Indeed,  for  small  taper  angles  it  could  be  useful 
« not  too  far  from  the  intersection  of  the  inner  and  outer  sur- 

faces as  it  stands.  But  it  would  be  better  to  try  to  develop 
the  idea  of  the  method  for  points  far  from  this  intersection. 

We  have,  indeed,  made  a stab  at  the  problem  and  hope  that  our 
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contributions  will  be  the  first  step  in  the  solution.  However, 
at  this  point,  we  feel  that  much  remains  to  be  done. 
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